FLAME PROPAGATION MODELING METHOD 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0001] The present invention relates to a flame propagation modeling method. 
Background Information 

[0002] One example of a flame propagation modeling method is described in an article 
"A Comparison ofFlamelet Models for Premixed Turbulent Combustion" that was 
presented in Combustion and Flame 95: 101 - 1 17 (1993). In this article, flame surface 
area density is defined as the flame surface per unit volume and serves as a method of 
modeling the flame propagation based on the fact that the flame progression transports, 
generates, and diffuses the flame surface area density. The flame propagation is estimated 
by modeling the flame generation in such a manner that the flame generation resulting 
from turbulent flow is inversely proportional to the chemical reaction characteristic time 
and proportional to the flame stretch rate. 

[0003] In view of the above, it will be apparent to those skilled in the art from this 
disclosure that there exists a need for an improved flame propagation modeling method. 
This invention addresses this need in the art as well as other needs, which will become 
apparent to those skilled in the art from this disclosure. 

SUMMARY OF THE INVENTION 
[0004] It has been discovered that with the aforementioned method, when the 
turbulence is small, the chemical reaction characteristic time becomes very large and the 
flame propagation is not reproduced. As a result, the flame propagation cannot be 
accurately predicted. 

[0005] Also, when the combustions of fields having different combustion modes are to 
be reproduced, it is necessary to multiply the flame generation corresponding to each 
combustion mode by a constant in order to adjust according to experimental values. 
[0006] In view of these issues, the object of the present invention is to provide a flame 
propagation modeling method that can accommodate a variety of combustion modes. 
[0007] Accordingly, in accordance with one aspect of the present invention, a method 
of modeling flame propagation is provided that comprises defining a flame surface area 
density of a flame as a flame surface area per unit volume of the flame; expressing flame 
progress as generation of the flame surface area density in terms of at least one of a 



turbulent combustion and a laminar combustion; determining flame growth resulting from 
turbulent combustion as being inversely proportional to a chemical reaction characteristic 
time and as a function of a turbulent Reynolds number; and modeling the flame 
propagation based on the flame growth. 

[0008] Accordingly, in accordance with one aspect of the present invention, a method 
of modeling flame propagation is provided that comprises defining a flame surface area 
density of a flame as a flame surface area per unit volume of the flame; expressing flame 
progress as generation of the flame surface area density in terms of at least one of a 
turbulent combustion and a laminar combustion; determining the flame growth resulting 
from laminar combustion as being proportional to both a laminar flame speed and to a 
ratio of a temperature of a burned portion to a temperature of an unburned portion and as a 
function of the Karlowitz number; and modeling the flame propagation based on the flame 
growth. 

[0009] These and other objects, features, aspects and advantages of the present 
invention will become apparent to those skilled in the art from the following detailed 
description, which, taken in conjunction with the annexed drawings, discloses a preferred 
embodiment of the present invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0010] Referring now to the attached drawings which form a part of this original 
disclosure: 

[001 1] Figure 1 is a simplified diagram that illustrates a basic concept of flame surface 
area density that is used in a flame propagation modeling method in accordance with the 
present invention; 

[0012] Figure 2 is a graph that illustrates a basic concept of flame growth during an 
initial stage of combustion that is used in a flame propagation modeling method in 
accordance with the present invention; 

[0013] Figure 3 is a graph that illustrates calculated values and experimental values of 

the pressure and heat generation in accordance with the present invention; 

[0014] Figure 4 is a diagram that illustrates a distribution of the flame surface area 

density; 



[0015] Figure 5 is a flowchart of an analysis performed when the combustions of 
fields having different combustion modes are reproduced in accordance with the present 
invention; and 

[0016] Figure 6 is a flowchart describing the setting of the analysis conditions in 
accordance with the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0017] Selected embodiments of the present invention will now be explained with 
reference to the drawings. It will be apparent to those skilled in the art from this 
disclosure that the following descriptions of the embodiments of the present invention are 
provided for illustration only and not for the purpose of limiting the invention as defined 
by the appended claims and their equivalents. 

[0018] Referring initially to Figure 1, a simplified illustration is shown to illustrate a 
basic concept of flame surface area density that is used in a flame propagation modeling 
method in accordance with the present invention. Since the actual flame is extremely thin, 
it is necessary to generate a calculation grid smaller than the flame band thickness in order 
to obtain a precise calculation. When such calculation grids are applied to an engine or a 
combustor, the number of calculation grids becomes enormous and the calculation cost 
becomes very expensive. By considering the flame surface area per unit volume, as 
shown in the Figure 1, it is possible to calculate the flame propagation precisely even 
when the calculation grid is larger than the flame band thickness. 

[0019] Referring to Figure 2, a graph illustrates the basic concept of the flame growth 
during the initial stage of combustion that is used in the flame propagation modeling 
method in accordance with the present invention. The horizontal axis indicates the 
amount of time that has elapsed since ignition (time after ignition) and the vertical axis 
indicates the flame radius. In Figure 2, the solid line St is the flame growth resulting from 
turbulent combustion and the broken line Sl is the flame growth resulting from laminar 
combustion. 

[0020] After combustion begins, the combustion first starts as a laminar flame. Then, 
turbulent combustion is gradually generated and eventually the turbulent combustion 
becomes large enough to ignore the laminar flame (e.g., when SI > S2). The growth of 
the flame is expressed by expressing the generation of the flame surface area density 27 in 
terms of the growth of the turbulent flame and the laminar flame. The generation of the 



flame surface area density E is expressed by combining the growth resulting from the two 
modes of combustion, i.e., the turbulent combustion flame growth St and the laminar 
combustion flame growth Sl. 

[0021] The turbulent combustion flame growth St is inversely proportional to the 
chemical reaction characteristic time and a function of the turbulent Reynolds number. 
Thus, the turbulent combustion flame growth St is expressed according to Equation 1 as 
follows: 

S T =a I (Ref 2 r^ Z (1) 
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where St represents flame growth resulting from turbulent combustion, 

E represents flame surface area density, k represents turbulence strength, 
e represents turbulence dissipation rate, Re t represents turbulent Reynolds 
number, r represents flame stretch rate, and ay and a 2 are model constants. 
[0022] The laminar combustion flame growth Sl is proportional to both the laminar 
flame speed Ul and to the ratio of the temperature of a burned portion to the temperature 
of an unburned portion (7yT w ) and is a function of the Karlowitz number. Thus, the 
laminar combustion flame growth S L is expressed according to Equation 2 as follows: 

S L = fa exp( - fcKa) £ U L Z 2 , (2) 

1 u 

where S L represents flame growth resulting from laminar combustion, 

E represents flame surface area density, U L represents laminar flame speed, 
T b represents burned gas temperature, T u represents unburned gas 
temperature, Ka represents Karlowitz number, and fij and f} 2 are model 
constants. 

[0023] Referring to Figure 3, a graph illustrates an example of an analysis using the 
flame propagation prediction model of the present invention. The horizontal axis indicates 
the crank angle after top dead center (degrees ATDC) and the vertical axis indicates the 
pressure (MPa) and the heat generation (dQ/d0). Both calculated values and experimental 
values are shown in Figure 3. Thus, it will be apparent to those skilled in the art from this 
disclosure that the flame propagation prediction model of the present invention can be 
used to optimize combustion in an internal combustion engine. 

[0024] Referring now to Figure 4, the distribution of the flame surface area density E 
is obtained from the example of the analysis of Figure 3 using the flame propagation 



prediction model of the present invention. The best result is obtained when the value of a 2 
in Equation 3 (as seen below) is 0.5 so that the square root of the turbulent Reynolds 
number is used. 

as + du£ = a Ujx) +ai(Re f2 r L £+ p iexp( .p 2Ka) ^ ULl 2_ D (3) 
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where E represents flame surface area density, k represents turbulence strength, 

e represents turbulence dissipation rate, Re t represents turbulent Reynolds 
number, /"represents flame stretch rate, Ul represents laminar flame speed, 
T b represents burned gas temperature, T u represents unbumed gas 
temperature, Ka represents Karlowitz number, v, represents turbulent 
kinematic viscosity, a c represents turbulent Schmidt number, 
D represents air resistance force, and a 2 , fii and f} 2 are model constants. 

[0025] Figure 5 shows a flowchart of the analysis performed when the combustions of 

fields having different combustion modes are reproduced. 

[0026] In step SI, the fluid analysis model (grid data) is acquired. In step S2, the 
experimental conditions (intake air temperature, air-fuel ratio, etc.) are acquired. 
In step S3, the fluid analysis grid is read in. In step S4, the analysis conditions are set 
according to the flowchart of Figure 6 (discussed later). 

[0027] With these steps, a calculation analysis grid is generated based on the analysis 
model and the ignition timing and initial analysis conditions which are set based on 
experimental values and other simulations. When the device being analyzed is an engine, 
the operation profiles of the intake valves and exhaust valves and the behavior profile of 
the pistons are set. 

[0028] The analysis condition setting flowchart of Figure 6 will now be described. 
The following processing operations are executed in steps SI 1 to SI 5: setting the intake 
gas composition, setting the boundary conditions, defining the physical properties, setting 
the initial conditions, and chemical reaction and combustion related processing. 
[0029] Returning to Figure 5, in step S5, based on the set initial values, the transport 
and diffusion of the gas flow and fuel inside the combustion chamber of the engine or 
combustor are calculated along with the fuel spray. The transport and diffusion of the gas 
flow and fuel are calculated by using such methods as calculus of finite differences, finite 
element analysis, and finite volume analysis and discretizing the fluid equations. The 
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turbulent flow is calculated using any of various turbulence models, such as direct 
numerical calculation, large eddy simulation, or a k-e model. In addition to solving the 
gas flow with good precision, the turbulence strength k, the turbulence dissipation rate a, 
the turbulent Reynolds number Ret, and the Karlowitz number Ka are calculated. 
[0030] The laminar flame speed Ul used is a value obtained experimentally or from a 
database or a value calculated using any of various laminar flame speed models developed 
by such persons as Keck and Gulder or a detailed chemical reaction model. 
[0031] In step S6, it is determined if combustion has begun, i.e., if the ignition timing 
set in step S4 has passed. If it is determined that combustion has begun, control proceeds 
to step S7. Meanwhile, if it is determined that combustion has not begun, control returns 
to step S5. 

[0032] In step S7, the transport of the flame surface area density L is calculated based 
on the aforementioned Equation 3 or the following Equation 4. 

where 27 represents flame surface area density, k represents turbulence strength, 

e represents turbulence dissipation rate, Re t represents turbulent Reynolds 
number, r represents flame stretch rate, U L represents laminar flame speed, 
T b represents burned gas temperature, T u represents unburned gas 
temperature, Ka represents Karlowitz number, v r represents turbulent 
kinematic viscosity, a c represents turbulent Schmidt number, 
D represents air resistance force. 
[0033] Similarly to the gas flow, Equation 3 or 4 is solved by using such methods as 
calculus of finite differences, finite element analysis, and finite volume analysis to 
discretizing the fluid equations. 

[0034] The flame stretch rate T is a coefficient expressing the stretch and quench of 
the flame and, here, it is determined using the ITNFS (intermittent turbulence net flame 
stretch) model. This model is a function of the ratio between the integral length scale 
developed by Poinsot et al and the flame band thickness and the ratio between the 
turbulence strength and the laminar combustion speed. It is a function that is formulated 
using a direct numerical simulation of the interaction between the flame and eddies, and 
experimental data related to intermittent turbulence and serves to provide the net flame 



spread. Also, the resistance force D that air exerts on the flame (flame resistance force) is 
found using Equation 5, which was developed by Poinsot et al. 
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U L +Ck l . 

D = y 1 2 (5) 

/Yfv 

where Y F represents fuel mass fraction, and k represents turbulence strength, and 
Yfv represents unburned fuel mass fraction. 
[0035] The size of the initial flame surface area density Z at the time of ignition is 
based on a predetermined flame kernel size and can be determined with Equation 6 using 
the fuel consumption mass niF during ignition. It is also acceptable to determine the flame 
kernel using an ignition model based on the air-fiiel ratio, turbulence strength, etc. 

CDFchem = m F U L £ (6) 

where co Fc hem represents fuel consumption rate, and m F represents fuel mass. 
[0036] In step S8, the ratio of the temperature of the burned portion to the temperature 
of the unburned portion (7VT U ) and the fuel consumption mass are found. Then, the field 
fuel consumption rate coFchem is found with Equation 6 using the transport of the flame 
surface area density S and the transport and diffusion of fiiel are discretized and solved 
using Equation 7 as follows. 

dY„ 



dt dx, dx: 



v, dY F 



V <T C dx, j 



\ 

+ COFchem + COpevap (7) 



where coFevap represents the increase resulting from vaporization of fuel mist. 
[0037] In step S9 after-processing is executed. The flame propagation inside the 
combustion chamber is predicted by executing steps S5 to S9 over a predetermined cycle 
time and number of steps. 

[0038] According to the present invention, the generation of the flame surface area 
density 27, which expresses the progress of the flame, is expressed with two modes, namely, 
the flame growth St resulting from turbulent combustion and the flame growth Sl resulting 
from laminar combustion, and the flame growth St resulting from turbulent combustion is 
inversely proportional to chemical reaction characteristic time and a function of the 
turbulent Reynolds number. Consequently, the turbulent flame resulting from turbulent 
combustion is strongly dominated by turbulent mixing and the combustion reaction 
proceeds in conjunction with the turbulent mixing. Since the turbulent mixing has a strong 
correlation with the turbulent Reynolds number, the flame propagation in various fields of 
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strong turbulence can be predicted by expressing the flame growth resulting from 
turbulent combustion as a function of the turbulent Reynolds number. 
[0039] Also according to the present invention, the generation of the flame surface 
area density Z", which expresses the progress of the flame, is expressed with two modes, 
namely, the flame growth St resulting from turbulent combustion and the flame growth S L 
resulting from laminar combustion, and the flame growth Sl resulting from laminar 
combustion is expressed based on the fact that it is proportional to both the laminar flame 
speed U L and the ratio of the temperature of a burned portion to the temperature of an 
unburned portion (Tf/Tu) and is a function of the Karlowitz number (exp(- fi 2 Ka)). 
Consequently, laminar combustion is assumed to be reasonably free of disturbances such 
that the combustion process is dominated by chemical reactions. Since the laminar 
combustion mode is characterized by the Karlowitz number, which is the ratio of the 
Kolmogorov time scale to the chemical reaction characteristic time, the flame propagation 
in various fields of weak turbulence can be predicted by expressing the flame growth as a 
function of the Karlowitz number. 

[0040] Also according to the present invention, the generation of the flame surface 
area density 27, which expresses the generation of the flame, is expressed as a combination 
of the turbulent combustion and the laminar combustion. Since the flame generation is 
expressed by combining the turbulent flame generation (which is expressed as a function 
of the turbulent Reynolds number) and the laminar flame generation (which is expressed 
as a function of the Karlowitz number), the flame propagation can be predicted even for 
combustion modes that are dominated by both laminar combustion and turbulent 
combustion. It is also possible to predict the flame propagation at fields where at the 
beginning of combustion the turbulence is very weak and laminar combustion dominates 
due to a very small Reynolds number, but eventually the turbulence generated by the 
combustion causes the turbulent Reynolds number to become large, strengthening the 
turbulent combustion. 

[0041] Also according to the present invention, the flame growth St resulting from the 
turbulent combustion is found based on the fact that it is inversely proportional to the 
chemical reaction characteristic time and proportional to both the turbulent Reynolds 
number raised to an exponential power (Re t a2 ) and the stretch rate r of the flame. 
Consequently, the turbulent flame is even more strongly dominated by turbulent mixing 



and the combustion reaction proceeds in conjunction with the turbulent mixing. Since the 
turbulent mixing has a strong correlation with the turbulent Reynolds number, a situation 
in which the flame speed increases as the turbulence of the field strengthens can be 
reproduced and the flame propagation in various fields of strong turbulence can be 
predicted by expressing the flame growth St resulting from turbulent combustion as the 
turbulent Reynolds number raised to an exponential power (Re t a2 ). 
[0042] Also according to the present invention, the transport of the flame surface area 
density, which expresses the generation of the flame, is expressed in two modes, namely, 
the flame growth S T resulting from turbulent combustion and flame growth S L resulting 
from laminar combustion, and the flame growth Sl resulting from laminar combustion is 
expressed based on the fact that it is proportional to the laminar flame speed U L , to the 
ratio of the temperature of a burned portion to the temperature of an unbumed portion 
(7VT W ), and to an exponential function of the Karlowitz number (exp (-yS 2 Ka)). 
Consequently, by expressing the generation of the laminar combustion as an exponential 
function of the Karlowitz number (exp(-/? 2 Ka)), it is possible to make the generation of 
laminar combustion large when the turbulence of the field is small and small when the 
turbulence of the field is large, it is possible to reproduce the flame growth S L resulting 
from laminar combustion and reproduce quenching, and the flame propagation in various 
fields of weak turbulence can be predicted. 

[0043] Also according to the present invention, the exponential function (exp(- fi 2 Ka)) 
is the base of the natural logarithm raised to the power of the Karlowitz number (- fi 2 Ka). 
Consequently, the flame growth Sl resulting from laminar combustion and quenching can 
be reproduced more precisely and the flame propagation in various fields of weak 
turbulence can be predicted. 

[0044] Also according to the present invention, the flame growth St resulting from 
turbulent combustion is expressed according to Equation 1 . Consequently, it is possible to 
reproduce a situation in which the flame speed increases as the turbulence of the field 
strengthens and the flame propagation in various fields of weak turbulence can be 
predicted. 

[0045] Also according to the present invention, the flame growth resulting from 
laminar combustion is expressed according to Equation 2. Consequently, it is possible to 
make the generation of laminar combustion small when the turbulence of the field is large, 



it is possible to reproduce the flame growth Sl resulting from laminar combustion and 
reproduce quenching, and the flame propagation in various fields of weak turbulence can 
be predicted. 

[0046] Also according to the present invention, the transport of the flame surface area 
density 27, which expresses the generation of the flame, is expressed with two modes, 
namely, the flame growth St resulting from turbulent combustion and the flame growth Sl 
resulting from laminar combustion, and the generation of the flame is suppressed by a 
resistance force D imposed by air. As a result, by expressing the generation of the flame 
such that it is suppressed by a resistance force D imposed by air, the effect of the air 
resistance D on the flame can be reproduced and the flame propagation can be predicted 
with good precision with respect to combustion modes that are dominated by both laminar 
combustion and turbulent combustion. 

[0047] Also according to the present invention, the transport, generation, and diffusion 
of the flame surface area density are expressed according to Equation 3. Consequently, 
the generation of turbulent combustion and laminar combustion and the resistance D 
imposed on the flame by air can be expressed and the flame propagation can be predicted 
with good accuracy with respect to combustion modes that are dominated by both laminar 
combustion and turbulent combustion. 

[0048] Moreover, terms that are expressed as "means-plus function" in the claims 
should include any structure that can be utilized to carry out the function of that part of the 
present invention. The terms of degree such as "substantially", "about" and 
"approximately" as used herein mean a reasonable amount of deviation of the modified 
term such that the end result is not significantly changed. For example, these terms can be 
construed as including a deviation of at least ± 5% of the modified term if this deviation 
would not negate the meaning of the word it modifies. 

[0049] This application claims priority to Japanese Patent Application No. 2003- 
066547. The entire disclosure of Japanese Patent Application No. 2003-066547 is hereby 
incorporated herein by reference. 

[0050] While only selected embodiments have been chosen to illustrate the present 
invention, it will be apparent to those skilled in the art from this disclosure that various 
changes and modifications can be made herein without departing from the scope of the 
invention as defined in the appended claims. Furthermore, the foregoing descriptions of 
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the embodiments according to the present invention are provided for illustration only, and 
not for the purpose of limiting the invention as defined by the appended claims and their 
equivalents. Thus, the scope of the invention is not limited to the disclosed embodiments. 



